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r-{ ■ In the hyperbolic slicing of de Sitter space appropriate for open universe models, a curvature 

scale is present and supercurvature fluctuations are possible. In some cases, the expansion of a scalar 
field in the Bunch-Davies vacuum includes supercurvature modes, as shown by Sasaki, Tanaka and 
Yamamoto. We express the normalizable vacuum supercurvature modes for a massless scalar field 
in terms of the basis modes for the spatially-flat slicing of de Sitter space. 
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Scalar fields in de Sitter spacetime have long provided a testing ground for issues of quantum field theory in curved 
Q" 1 spacetime Further motivation for their study stems from the central role they play in inflationary cosmology 

|3j. Several different coordinate systems can be used to cover de Sitter space, and subtleties in the quantization of 
fields can arise in some of the less familiar coordinate systems. These subtleties have been highlighted by recent 
models of open inflation ^J^], in which two periods of inflation are separated by nucleation of a bubble. The bubble 
interior includes an open universe (Qo < 1), where we could be living today, described by hyperbolic, spatially-curved 
coordinates ||. 

A key difference between these hyperbolic coordinates and the more familiar spatially-flat slicing of de Sitter space 
is the presence of a curvature scale. This in turn leads to the possibility of supercurvature Q fluctuations, fluctuations 
with wavelength longer than the curvature scale. Unlike the continuum of modes familiar from the spatially-flat slicing 
of de Sitter space, a normalizable supercurvature mode may exist for an isolated, discrete eigenvalue of the spatial 
Laplacian, or not at all. Although such modes have no analogue in the spatially-flat slicing of de Sitter space, it 
has been shown Q (see also |j) ) that the supercurvature modes must be included in the vacuum spectra of low-mass 
scalar fields in order to produce a complete set of states, and hence the proper Wightman function. 

For massless, minimally-coupled scalar fields in de Sitter space, there is in addition a well-known infrared divergence 
in the (coordinate-independent) Wightman function. The infrared divergence is related to a dynamical zero mode in 
the spectrum of a massless field • Kirsten and Garriga [|ll| covariantly quantized this zero mode in a spatially- 

closed slicing of de Sitter space. Extending their result from these closed coordinates to the coordinate system 
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appropriate to open inflation has not yet been done. Here we identify the zero mode as one of the supercurvature 
modes when quantizing the minimally-coupled massless field in the open hyperbolic coordinates. 

It has already been noted that supercurvature modes can make significant contributions to the fluctuations in the 
cosmic microwave background (CMB) radiation, and several of their effects have been calculated for models of open 
inflation |l2| (The massless field zero-mode subtleties, except for a variant studied in |L5|| , do not arise for these 
specific CMB calculations, which are sensitive to higher multipoles.) In addition to contributing to observable density 
fluctuations, such long wavelength, supercurvature modes might play a role |]l6f at the end of inflation, when recent 
advances |fl7| in the theory of reheating are taken into consideration. 

In summary, increased understanding of these supercurvature modes is motivated both by general questions of 
quantizing fields in curved backgrounds, and by recent inflationary model-building. We will focus here on the example 
of supercurvature modes for a massless, minimally-coupled scalar field, expressing them as a sum over the basis modes 
for a spatially-flat slicing of de Sitter space. This overlap gives a measure of "where all the supercurvature modes go" 
in the familiar flat-space spectrum of such fields. The massless case is chosen for tractability, and questions about the 
zero mode are postponed for future work Jl8{ . 

Throughout this paper, we consider only an unperturbed de Sitter metric; a more complete treatment would include 
study of the backreaction of such fields on the background metric. Because the supercurvature modes stretch beyond 
the horizon, any such study of the coupled metric fluctuations would need to pay special attention to the gauge 
subtleties which always accompany superhorizon fluctuations fl9| ], and such issues are not pursued here. 

In section II, the two covers of de Sitter space (including the flat and hyperbolic slicings) are given, and the 
field quantization pertinent to open inflation in both systems is reviewed. As these supercurvature modes have no 
analogue in the usual flat slicing of de Sitter space, this section gathers some previous work on supercurvature modes 
and provides notation and context for the rest of the paper. Section III specializes to the massless case. The explicit 
calculation of the overlap for some of these supercurvature modes and the more familiar flat space modes is given, 
and indicates a general form for the overlap between all the massless supercurvature and spatially flat modes. We 
verify this general form by integrating over the spatially flat modes, weighted by the overlap, to obtain the original 
supercurvature modes. Concluding remarks follow in Section IV. Three Appendices include the supercurvature mode 
normalization at fixed time in the flat coordinates, and details of the overlap calculation along different hypersurfaces 
within de Sitter space. 

II. DE SITTER SPACETIME AND SCALAR FIELD QUANTIZATION 

We begin by considering de Sitter spacetime as embedded within a 5-dimensional Minkowski spacetime, with the 
five coordinates subject to the constraint pj: 
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-(-°) 2 + X>0 2 = i- a) 

2 = 1 

The radius of the embedded spacetime, H^ 1 , is scaled to unity. 

A spatially flat slicing (see, e.g. which partially covers the resulting 4-dimensional de Sitter spacetime is 

z° = Binht, + \e^r) = ~ (l + r}) + ~rj, 
z 4 = cosht/ - -e^r} = ~^~ rj) - \v, 

F=e^r/fi = -^0, (2) 

where Q, = (sin#cos</?, sinfsin^, cos(9), and conformal time drj = a~ 1 dtf 1 or ri = — a (tf) — —e~ tf . Here r/ > 0, 
— oo < rj < 0. (See Figure 1.) 




FIG. 1. De Sitter spacetime as covered by the coordinates in equation (g). Solid lines are lines of constant r, and dashed 
lines are lines of constant r\. 

The metric on this portion of the spacetime is 

ds) = rf 2 [-dri 2 + drj + rjdn 2 ] , (3) 

where df2 2 = d0 2 + sin 2 9d(p 2 . The coordinates xj and metric ds 2 are useful in ordinary models of inflation, in which 
the spatial curvature quickly becomes completely negligible. These coordinates xj cover only that half of the total 
spacetime with z° + z A > 0. Replacing z 4 — * — z A in the above produces a second flat coordinate system. 

In contrast, for models of open inflation, open hyperbolic coordinates are appropriate inside the open universe. 
This coordinate patch is part of the full de Sitter spacetime (discussed in detail by [R^Of) as shown in Figure 2. 
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FIG. 2. De Sitter spacetime as covered by the coordinates in equation (|J). Solid lines are lines of constant r, and dashed 
lines are lines of constant t. If a (thin-walled) bubble nucleated with center at point P, then region R, the forward light cone 
of point P, would contain in the interior of the nucleated bubble. In models of open inflation, region R would contain our 
observable universe at some time t inside the light cone corresponding to t — 0. 

The two most important of these regions for our purposes here are region C, a large, compact subspace, and region 
R, an open hyperboloid bordering region C. These are related to the embedding coordinates by 

cos t c sinh r c = sinh t r cosh r r , 
sint c = cosh£ r , 

cos t c cosh r c f2 = sinh t r sinh r r fl. (4) 

7T 

— — < t c < — , — oo < r c < oo and < t r , < ry, (5) 
and are related by the analytic continuation t r = it c — ^ and r> = r c + ^ ■ The corresponding metrics are 

ds 2 c = dt\ + cos 2 t c [-dr 2 c + cosh 2 r c dn 2 ] , 

ds 2 r = -dt 2 . + sinh 2 t r [dr 2 r + sinh 2 r r dVL 2 ] . (6) 

Note that r c plays the role of 'time' within region C. 

The forward light cone of the point P in Figure 2 is the center of the nucleated bubble in models of open inflation, 
so that region R contains the spatially-open universe we may be in today. Region C is of interest because surfaces of 
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These coordinates lie in the ranges 



4 



fixed time, r c = constant, correspond to Cauchy surfaces for de Sitter space. [J] Quantizing on a Cauchy surface in C 
thus specifies all of the "initial data" for the system, including the initial conditions for the spatially-open universe 
inside region R. 

The equation of motion for a scalar field of mass M in this space is 

(□ - M 2 ) 0(a) = ((^r 1 ^ [W^gTdv] - M 2 ) fa) = 0, (7) 

and separation of variables gives a family of solutions 4>kim(x). To quantize, (f> and its canonically conjugate momentum 
n are promoted to Heisenberg operators, and expanded as 



»> = E[* 



| 4>kt m (x)a H m + <t>ktm( x )&lim ( 8 ) 

k£m 

and similarly for EL The creation and annihilation operators satisfy [autm, ^k'l'm'] = $(k — k')8in8 m 'm with the other 
commutators vanishing. Here an overline denotes complex conjugation, ^2 kem is a placeholder for the appropriate 
measure, and the choice of vacuum state satisfying ftfe^ m |0) = for all (k,£,m) provides a division into positive and 
negative frequency modes. The fixed-time canonical commutation relations, [<j>(x), Ti{x')] = iS 3 (x — x') then imply 
that 

(fom, <t>k>l>m') =S(k — k') 5e>i5 m > m , {4>Mm> <t>Wt'm') = ' 

(fyklm, <t>k't'm') =S(k- k') St'iSm'm, (9) 

where the Klein- Gordon inner product is defined by 



(4>klm, 4>k'£'m') = ~ J \4>ktm ^ $ k'i'm' J V~SS dY. . (10) 

Here E is a (spacelike) Cauchy surface and n M is a future-directed unit vector normal to this surface. (The extra factor 
of i which multiplies the righthand side of equation ( |io|) in jj| is absent here because of our different sign convention 
for the metric.) This inner product is independent of Cauchy surface, and more generally is independent of the choice 
of E, as long as the fields fall off sufficiently quickly on the time-like boundaries. 

The physically-motivated choice of initial vacuum state in models of inflation is the Bunch-Davies vacuum [p2[ , 
which respects the symmetries of de Sitter space and reduces to the Minkowski space vacuum at early times and over 
short distances. For the flat slicing, the Bunch-Davies positive-frequency modes are 

fornix,) = -$=exp + ^^H^{-kri)u{kr,)Y lm {n), (11) 
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1 A Cauchy surface is any hypersurface such that every future-directed time- like vector intercepts it exactly once. (See, e.g. 
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where v=J\- M 2 , hP(z) is a Hankel function of the first kind, je(z) is a spherical Bessel function, and Ye m (il) 
is the usual spherical harmonic. The measure in equation (|^) for the expansion in flat modes tfikim is Ylkem = 

Jo aK l^i=0 L^m=-t 

For the hyperbolic slicing with the metric ds 2 of equation (^|), the positive- frequency solutions to the equations of 
motion are p|,p|,p3| 

u P tm{x) = *Y^f P i{r c )Y lm (n), (12) 

<X[t c J 



where 



P ^ i Vn 2 (i sinhrv) 

f pe (r c ) = IP ~ 12 \ \ (13) 
\Ji cosn r c 



and 



v' = V-\ = \hr-M 2 -l . (14) 
2 V 4 2 v ' 

Here P^(z) is an associated Legendre function of the first kind, while the specific forms of a p ,(3 p will not be needed. 
Within region C, the f p i(r c ) play the role of positive-frequency solutions, and the Xp(^c) are spatial eigenfunctions; 
when continued into region R, these roles are reversed. 

A Friedmann- Robertson- Walker metric with spatial curvature K and cosmic scale factor a(t) has a physical curva- 
ture length scale a{t)/\K\. For a flat universe the comoving curvature length scale thus runs off to infinity, whereas 
in a spatially closed or open metric, the comoving curvature length scale is +1. Eigenvalues of the (region R) spatial 
Laplacian in this background are — (fc/a) 2 , where (k/a) is the inverse of a physical length: k/a — l/xphys, with 
< k 2 < oo. Defining p 2 = k 2 — 1, p 2 > for subcurvature modes, and —1 < p 2 < for supercurvature modes. 
The continuum of subcurvature modes < p < oo is sufficient to describe a Gaussian random field in region R, 
see for detailed discussion. In addition, for p 2 < 0, inner products of the form (|l(]) on fixed-time (non-Cauchy) 
surfaces within region R diverge, and so all supercurvature modes naively appear to be unnormalizable. Studying 
quantization and completeness more appropriately on a fixed time r c Cauchy surface in region C, it was found in 
H that for M 2 < 2H 2 (restoring the Hubble radius, H~ 1 ), supercurvature modes are normalizable in vacuum for 
a discrete value of p. In addition, it was shown there that this value of p must be included to obtain the correct 
Wightman function for the Bunch-Davies vacuum. (See also Q.) 

This discrete normalizable supercurvature mode can be understood as follows |^,^3|. Its presence is suggested by 
an analogy between these supercurvature modes and bound states in a potential. The only dependence on M in the 
wavefunctions is in the 'spatial' (in region C) eigenfunctions Xp(^c)- Defining smt c = tanhw, the equation of motion 
for Xp becomes 



G 



Xp 
U(u) 



p 2 x P 



M 2 — 2 
cosh 2 u 



(15) 



This is a one-dimensional Schrodinger-like equation with the potential U(u) and energy p 2 . As noted in |lqjr4| , p3| , 
the potential U (u) vanishes as u — > ±00, revealing that in this limit there exists a continuous spectrum of modes with 
p 2 > 0. But over finite intervals of u, if the mass of the field satisfies M 2 < 2, the potential U(u) has a valley and 
the modes \p behave as discrete bound states with p 2 < 0. As DeWitt has shown, such discrete states in a field's 
spectrum are generic for fields quantized on compact subspaces. p3 
The inner product on this space is proportional to 



du x p (tanh u)x p i (tanh u) 



(16) 



so for normalizability, x p (tanhw) must be bounded as u — > ±00. (A similar argument is found in the appendix of |p[.) 
A supercurvature mode has p 2 < and p 2 + 1 > 0, so define p — iA, with < A < 1, and A real. The asymptotics of 
X P near tanhu — > ±1 yields (cf. p^]): 



r(l T tp)P^ lp {t&nhu) ~ e ±lpu , 



(17) 



This is finite only for PJi p since ip = — A < 0, thus (3 P = for the supercurvature modes. The limit of P,/ (tanh u) 
asm -00 is |2lR 



T(l + A)P^ A (tanh 



r(i + A)r(A) A|M| r(i + A)r(-A) 



e A l"l , u 



-00. 



(18) 



r(l + v' + A)T(A - v') T{-v')T{l + v') 

As < A < 1 is non-negative and < Re v' < 1, the coefficient of e A l"l vanishes only if A = v' ', producing an isolated 
value of p which is normalizable. (Note that for subcurvature modes, with p real and non- negative, the solution Xp(^c) 
simply oscillates at both endpoints and so remains finite.) 

The Klein Gordon normalized supercurvature modes are then 

-£-1/2 



uum{x c ) = N Ae . m (-icost c ) 

NAim 



p 

-1 -A-l/2 



(i sinh r c ) 



V^cosnr c 

T(v' + l/2)T{-u' + £+ l)r(i/ +£ + l)" 1 1/2 



20FT(i/) 



(19) 



The scalar field in region C is expanded as 



2 Note that in the derivation of this result in their equation (2.9) is incorrect, including only the factor w M instead of 
[(1 — w)/w] 11 , though their next equation is correct. 
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poo 

4>(x c ) = I dp^2 [upi m {x c )a p t m + H.c] + ^2 {uMm{x e )aum + H.c] 

J° i,m t,m 

= $(p)(x c ) + <P(A){x c ), (20) 

where "H.c" denotes the Hermitian conjugate, and Yltm = ^2iLol2m=-e- Once quantized, these modes may be 
continued into region R. In the presence of a bubble wall rather than in the vacuum, the value of A may change 
and a supercurvature mode may appear even for M 2 > 2, depending on the details of the model. The normalizability 
conditions may be solved for numerically. In the absence of gravity, there is also a supercurvature mode for the 
fluctuations of the bubble wall itself, which appears to become singular once gravity is included |26| ] . 

III. OVERLAP OF SUPERCURVATURE WITH FLAT MODES FOR A MASSLESS FIELD 

We now re-express these supercurvature modes UA=ifon in terms of the spatially-flat modes <pktm for a massless 
scalar field. There are two things to consider, however, before proceeding to the calculation. First, for M = 0, the 
£ = supercurvature mode of equation ( |l9|) is constant. Consequently its Klein-Gordon inner product is zero and 
its norm (proportional to the square root of the inverse of this norm) diverges as Y x l 2 (£ = 0). Including this state 
naively in the Wightman function sum over states will thus diverge as well. It is known that the Wightman function 
for a massless, minimally- coupled scalar field in de Sitter spacetime is infrared divergent |To|Jri|| , and in this particular 
slicing, the specific state ma=i,oo appears to be the lone source of the divergence. 

The r c -dependent portion of the equations of motion in region C, are exactly like a (2, l)-dimensional spacetime 
with scale factor a(r c ) = coshr c , For positive frequency modes, 



1 9 ( 2, , d \ , . , ~ , £(£+1) 



a 2 (r c )dr c \ dr c J a 2 (r c ) 

and the 'frequency' uj Tc associated with the time-coordinate may then be written 



{f P i(r c )Ye m (n)} = (21) 



< = (l+p 2 ) + ^±^, (22) 



the analogue of u 2 = M 2 + k 2 /a 2 (t). This Ctv e vanishes only for p = i and £ = 0, indicating a zero mode corresponding 
to the symmetry <j> — ► (f> + constant. This zero mode should be replaced by a collective coordinate.^ [ pl] , ^8|j27| Zero 
modes have also been identified in the context of two field models in . In the following, we address only the finite 
modes (£ > 0) for the massless, minimally-coupled case. 

The only Cauchy surface for the entire de Sitter spacetime in region F is the limiting curve 77 — > 0, corresponding to 
t r — > 00. This is in contrast with region C where any 'time' r c — constant surface is a Cauchy surface. However, all of 



3 Such divergences do not directly affect CMB anisotropy calculations, which correspond to £ > 0. 
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region R, containing our open observable universe in models of open inflation, is contained within region F. Because 
our aim is to provide a useful heuristic relation between the unusual supercurvature modes and the more familiar 
spatially-flat modes, we will thus work in this section at fixed time within region F. Verifying these results, Appendix 
B contains a parallel calculation for all odd £ and I — 2 along the proper Cauchy surface r c = 0. In Appendix C the 
overlap for £ = 1 is found on the boundary of F, corresponding to Cauchy surface 77 = 0. 
The Klein-Gordon inner product in region F is 



(u, v) — —ia (77) J drfTf J dCl [u (d v v) — (d n u) u] s . (23) 

This can be evaluated along any fixed-77 surface £, and, if the fields fall off sufficiently quickly on the time-like 
boundaries, will be independent of the specific choice of 77, even though such fixed-77 surfaces are not Cauchy surfaces 
for the entire spacetime. The modes 4>kem(xf) and <fikem( x f) satisfy the inner product relations of equation (^J) along 
such fixed-r? surfaces within region F, and form a complete set of orthonormalized modes within region F; thus they 
may be used to expand any normalizable function within region F. The massless supercurvature mode «A=i/'m' is 
normalizable on fixed 77 surfaces with — 1 < 77 < 0, as shown in Appendix A. Its norm being independent of 77 suggests 
that its falloff is fast enough to make the inner products independent of 77 as well. Thus we expect we can express the 
massless supercurvature mode u^ = \^ m i (x r (xf)) in terms of the flat basis functions as 



i,m 

with oiktm and Pkim constant complex coefficients. 



/■oo 

(x r (x f ))= I rffoV] \ctklm 4>klm(xf) + Pkfrn 4>Mm( X f)] > ( 24 ) 

Jo ™ 



OLkl-n 



— (UA=l,em, <t>klm) , fiklm — — {uK=ll m , 4>ke m ) ' (25) 



even for a fixed 77 (non-Cauchy) surface. 

For M = 0, the spatially-flat basis modes (equation (111])) reduce to 



cj>Mm(xf) = —=(l + ikr])e lkv j e (kr f )Y em (n) 
V nk 

F ke ( V ,r f )Y em {Q), (26) 



and the normalized supercurvature modes for the massless case are 

p ^ 1/2 / cosh v \ 

UA=Mm(av) = -^r(l)r(l + 2) ~ 3/2 ^ —Y tm {Q) 

£ V sum TV 

1 



= -^/T(£)T(e + 2)Sfir r ) Y £m (Q) (27) 

within region R. Using the embedding coordinates of the five-dimensional Minkowski space to relate the coordinates 
Xf and x r (see equations (||) and (|J)) yields 
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where 



The coordinate g is convenient because of the identities (see |2^], equations 8.2.7 and 8.6.7): 



(29) 



P 



1 



)\{z) = V^l 



r(a + /3 + l) 
dQt(z) 
dz 



(30) 



where Re [z] > 0. Taking g = z and noting that < r/ < oo in the inner product means that we want fixed 77 surfaces 
with — 1 < T] < 0. When z 2 = g 2 < 1, the argument of PZp—in becomes complex and hence the righthand side 
should be understood with z having a small imaginary part. Unless r\ = 0, fixing 7/ and letting r/ range over its values 
< Tf < 00 will include some values of g 2 < 1. 
With the identities (BOh above, 



-P_3 /2 1/2 (cosh r>) 



V / smhrv 

and the overlap in equation (pq) is then 



1 



7rr(£ + 2) 



(, 9 2 - 1) d g Q e (g), 



(31) 



a Mm = -Ckii-rj) V fc " / dr, r 2 f j e (kr f ) [k 2 rj S e (r r (x /)) - (1 - iAt7)ft,ft(r r (x/))] : 



r(£)r(£ + 2) 



7rA: 



1/2 



(32) 



Note that the constant coefficients are purely real. Similarly, 

/3 Mm = C M (-r/)- 2 e- lfef '(-l) m 5 m ,_ m - [fc 2 ?7 - (1 + ik V )d v ] / ^ [r 2 j e (kr f )S e (r r )] (33) 

Jo 

As both <fi(xf),u(x r ) correspond to positive frequency for the Bunch-Davies vacuum, Pkim — for all k, £, and m, 
giving 

k 2 r\ 



/•OO 

^( fc )|s = / *7 [ r / ji( kr f)Si(r r )] s , 



(34) 



-2iC fe ^ 3 



a ikr/ 



{1 + ikr)) 



h{k)\ 



(35) 
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Equation ( p4| ) may be used to demonstrate explicitly that d v ctkim — identically allowing a choice of any convenient 
value of 77 to evaluate h(k). We choose r] = — 1 in this section; the case 77 = and £ = 1 is in Appendix C. 
Along the surface 7/ = —1, = r//2 and 

«^ " -^TWW f « - !) J ' (tt " ) ^' (36 > 

where the * indicates that a particular value of 77 has been chosen for this evaluation. For several small values of £, 
repeated integration by parts gives the general form (dropping the subscript * ) 

J dg g 2 (g 2 - l) j e (2kg)d g Q e {g) = a x cos(2fcg) + a 2 sin(2fcg) 

+ a 3 (ci [2fc(l + <7)]-ci [2fc(l-0)]) 
+ a A (Si [2fc(l+ 9 )]-Si [2k(l-g)]) 

+ a 5 Si[2kg], (37) 
where ci (z) = — dt £2IlL i Si (z) = J Q dt S! f^-, and generally at — ai(g, fc). For I = 1, the nonzero coefficients are 

ai 
a 2 



5F + 8FH^Tj [1 + t2 - MV1 ' 

(jVfj> + ' Jill) [2, + 

4fc 4 8fc 4 \g-lj 1 y n 

° 3 = sis t( fc2 ~ ^ cos ( 2fc ) ~ 2fc sin(2fc)] , 



8fc 5 
1 



0.4 



[2fe cos(2fc) + (fc 2 - 1) sin(2fc)] . (38) 



8fc 5 

Using |^5|,^8), as g — > 00, ax — > and 02 — > (6fc 2 ) -1 and thus 

/"°° 1 7rd — ik)e ik 

J dg g 2 (g 2 - l) j 1 (2kg)d g Q 1 (g) = ^ sin(2fc 5 ) ff ^ 00 + V 4fc5 7 [fc cos(fc) - sin(fc)] . (39) 

Based on comparison with the explicit calculation of this definite integral along different surfaces (see Appendix B), 
we drop the first term, as its limiting value oscillates at g — * 00. Then the coefficient of expansion Qtk,£=i,m becomes 

a M= i, m = -2ifc- 1 /2j 1 ( fc ). (40) 
Repeating the same analysis for £ = 2,3, for which 01 and a 2 both vanish identically as g — > 00, yields 

ak,e=2,m = -2iV3 k~ 1/2 j 2 (k), 

«M=3,m = fc- 1/2 j 3 (fc). (41) 

These first three expansion coefficients are plotted in Figure 3. 
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i overlap (k, ell) 



ell = 1,2,3 




FIG. 3. The first three expansion coefficients, iatem for I = 1, 2, 3, plotted against wavenumber k of the modes in the 
spatially-flat slicing, xj. 



(42) 
(43) 



These coefficients are finite in both the k — > and k — > oo limits, 

lim fc _1/2 j £ (fc) - constant k e ~ 1/2 -> (for ^ > 0) 

lim AT 1/2 jf(fc) — constant fc~ 3 ^ 2 sin ( fc - — J — > 

by the asymptotics of spherical Bessel functions {e.g. pq] , equations 10.1.4, and 9.2.5). 

Equations (^,^l|) suggest that the general form for the overlap of the supercurvature modes with the fiat ba- 
sis functions is proportional to je(k)k^ 1 ^ 2 . This can be tested by seeing if these postulated otkim reconstruct the 
supercurvature mode, i.e. 



UA=i,e m (x r ) 



dk aum 4>Un 



(44) 



The most convenient choice of r\ for this integral is the surface r\ — 0, since in this case the k dependence in <pe m is 
proportional to jt{kr f)k~ x l 2 . Writing aktm = T^e fc" 1 ^ 2 ji(k) and substituting into equation ((44|), we have 



-^T(£)T(e + 2)S e (r r (v = 0)) = ? -= V e 

2 \ 7T 



dk: 



je(kr f )j e (k) 



For the left hand side of this equation, using equation (|29| ) and setting r\ = 0, 



S e (r r (rj = 0)) = 



' 1 — r 



2r, 



-3/2 



1 + r / 



(45) 



(46) 



Region R corresponds to z 4 > 1, which requires r 2 < 1 as 77 — > 0. 
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The right hand side of equation (E^) can be integrated to give ( ]25| ], equation 11.4.34): 

, fe M^^ = ^_|iW. F(< ,_ 1/2;t + 3 /2;r? ), (47, 



(I 



where F(a, b; c; z) is the hypergeometric function. By using F(a, b; c; z) = F(b, a; c; z), the representation of Py{z) in 
terms of hypergeometric functions ( p8[ , equation 8.772.3), and Pj^{z) = P l ^ v _ l (z), the righthand side of equation 
(fl5|) becomes 



^ r[£W r{£ /2) p -/ -1/2 / 1+^ 

4 T^ + 3/2) V 7 ; f V / - 3 / 2 \l-rj 



'7T, 



' ~ r / D-^-l/2 [ 1 + r / 



V< f »V^T^ l^fj' (48) 

Thus the r/ dependence on both sides of equation ( |45| ) matches exactly, and the constant 2?£ may be read off: 

V e = -iy/2£(£+l). (49) 

This coefhcient reproduces the specific akim calculated earlier for i = 1, 2, 3 along the surface rj = — 1. Resumming 
the equation for r/ > 1 would correspond to region L in de Sitter space. 

Thus we conclude that the constant coefficients akim which relate the normalized supercurvature modes u\im and 
the spatially-flat basis modes <f>kim are: 



a kim = -iy/2e(t+l) k- 1 ' 2 Je (k). (50) 

This is the main result of this paper. 

Taking the limit r/ — > 1 in equation (^7|), and using the identity ( Ipgfl , equation 9.122.1): 

P , fl n r( 7 )r( 7 -q-/3) 

F(a, 0, 7; 1 = — — -jr , (51) 

r(7 - a)T(7 - (3) 



it is easy to verify in addition that 



,lk\n, l „ l \" = 2m+ li / rffc = 1, (52) 



JO ^ 

IV. CONCLUSION 

In conclusion, we have given the explicit form for the overlap between the flat basis functions (p^ ) and the massless 
supercurvature modes. As a result, the supercurvature modes within that patch of de Sitter space which would contain 
our open observable universe can be written 
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UK=H m {x r {Xf)) = —iy/2£(£+l) / dk k 1/2 j e (k) (j>Um{xf). (53) 





The long-wavelength supercurvature modes are distributed over the spatially-flat basis modes, oscillating over flat- 
space comoving wavenumber k with decreasing amplitude. More quantitatively, the spherical bessel functions je(k) 
have their first and largest maximum ( |2J| 10.1.59) near k ~ (£+ |)[1 + 0(£~ 2 / 3 )] with the approximation improving 
as I increases, but the damping envelope going as fc" 1 / 2 lowers this peak for higher I. It may be possible to use the 
description jl5| of M 2 ^ supercurvature modes as small perturbations of the massless supercurvature modes to 
extend the above to small M 2 . 

This expression for the supercurvature mode on fixed r\ surfaces, extending into region R, may be useful for 
understanding the effects of supercurvature modes during reheating. Unlike the event of bubble nucleation, reheating 
occurs in the future of region C and hence descriptions using fixed time r c surfaces are not appropriate. Rather, 
it is important to understand the dynamics of these modes within region R, corresponding to our observable, open 
universe. Having an expression for the normalizable supercurvature modes on slicings extending into region R is a 
step in separating long wavelength properties of these modes from issues related to their non-normalizability at fixed 
time t r in region R. 

We showed by comparing Cauchy and non-Cauchy surface calculations that in some cases non-Cauchy surface 
calculations of norms (in the appendix) and overlaps (in the text) agree for supercurvature modes, up to an identifiable 
boundary term. These non-Cauchy surfaces (fixed time in the flat coordinates) extend into the open universe and thus 
could be used (with caution) to calculate other properties. The complementary calculations along different surfaces 
were required to verify that the non-Cauchy surface representation was indeed correct. 

In addition, we identify one specific supercurvature mode as responsible for the well-known infrared divergence for 
massless scalar fields in de Sitter space. Not only does this mode have divergent norm; it is demonstrated here to be a 
dynamical zero-mode. This identification is a necessary first step toward its eventual replacement with an appropriate 
collective coordinate, similar to what has been done in closed coordinates JO]. Something similar has been done in 
fl5f in the context of two field models and quasi-open inflation, here it is found more generally as a property of the 
massless supercurvature modes. 

V. APPENDICES 
A. Normalization of supercurvature mode at fixed time — 1 < rj < 

In region F the supercurvature mode can be written as (using equation |28| ) 

cr t\\ -P-3/2 1 (coshr r ) D -l-i/2 ( 9 \ 

Se{rr{9)) = Vsinh, r = {S - 1} P -3/3 I 77=1 J ( } 
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which can be used to extend out of region R, to where g < 1. We will consider only — 1 < rj < for convenience. 
The inner product for fixed time r\ is 



oo 



(S e , S e ) = -ia 2 {rj) J dr f r 2 f J cM [S e (d v S e ) ~ (d v S e ) S t ] E . (55) 

The integral over 17 gives a delta function and will be suppressed. 

For g > 1, Si = Si because both the associated Legendre function and its argument arc real, and so the integrand 
disappears. For fixed r/, g 2 = 1 corresponds to 

(1 - r f f = => r f = 1±7] . (56) 

With — 1 < T] < 0, the integral is 

(S e , St) = i- [ dr f r f [S e (dgS t ) ~ (d g Se)S e ] s . (57) 
V Ji+n 

where we have substituted as well d v f(g) = —(rj/r f)d g f(g) and a 2 (?y) = rf 2 . In the region of integration, we also 



have Se(g) — Si(—g). Inside the integral, the terms where the derivatives act on (g 2 — l) 1 / 4 and its complex conjugate 



cancel out. Defining y — gj g 2 — 1, 



S t (d v Si) - (d n Se)S e = \g 2 - 1|^ f [PZ^iv) (d v Plfc'\-yj) 



9 y P:t 2 1/2 (y)) Plt 2 1/2 (-y) (58) 



dg 

':>■' .'2 / (yj) ' .; 

So we are left with the wronskian of PZ3J2 iv) an d PZ^J^^i - V) ti mes dy{g)/dg. We can now use ( equations 
8.736.2, 8.334.3, 8.335.1, and ||, equation 8.1.8) 

P^^i-v) - ^(y)P^(-y) = + M^[P^v)-^Qi{v) - ^(y)Q^y)] 

2 2 2 AT( M±^±2 )r( ££±^±l) 

= sin[(i/ + fi)7T\ 



7T (1 -y 2 )T{ ~^ +2 )r( ) 

2 2i TT 1 1 ' J (l-y2) r (^±l)r(|) 

2 sm[-(^ + 2)7r] r(_ ^ _ " 



(1-J/ 2 )7T L v T(£) 

2 1 



(l-y 2 )T(£)T(£ + 2) 



(59) 



Including 



£v(s) = ~ {g2 _! 1)3/2 > (6°) 



using |g 2 — l| 1 ^ 2 /(g 2 — l) 1 / 2 = i, and substituting the above, the inner product becomes 
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n jj v drf rf [Se ^ - {daSi) ^ = = - T{i)T{t + 2) \ LI drf rf 



1 



(l-r?) 2 + 



2 1 



4 

(61) 



r(£)r(^ + 2) 

which is independent of rj as promised. This also shows that the supercurvature modes are properly normalized 
for fixed time surfaces — 1 < r\ < in region F. This suggests that at spatial infinity in the flat coordinates the 
supercurvature modes have sufficiently fast falloff to correspond to their inner product on a Cauchy surface. (This is 
not true for fixed time surfaces in region R for example, as the supercurvature modes diverge at spatial infinity.) 

B. Overlap along r c = 0, for I odd and £ = 2 

In region C the inner product (equation (|l0|)) is 



dt c cosi c / dfl I — — v — u— — I (62) 

and will be taken on the Cauchy surface r c = 0. From their definitions, equations (||) and (Q), the spatially-flat 
coordinates and the region C coordinates are related via 

cos t r cosh r r 



'7 



sin t c + cos t c sinh r c 



V = ~ . , , 1 + • , , (63) 
sin tr + cos t r sinn r r 



with t c > 0. For 7j c < we need the analytic continuation of the basis function 



<t>Um{xf{r c ,t c )) = ^ (I + ikv(r c ,t c )) e-* k ^^Mkr f (r c ,t c ))Y e Un)- (64) 

V 7TK 



As this basis function has no branch cuts as a function of rj(r c , t c ), rj(r c , t c ), the continuation to t c < is straightfor- 
ward and the same for both positive and negative frequency mode functions. 

Using 4>kim{x) to denote both the function in region F and its analytic continuation into region t c < 0, the inner 
product ctkim between ma=i em and 4> ls 

,7r/2 , 

{uh=Hm,(l>U'm>) = —i / dt c COSt c / d£l 



'-tt/2 

x [u em (r c ,Q) (d r J> ulm ,(ri,rf,n)) - {d rc u lm {r c ,Q))4> M , m ,] \ Tc=0 . (65) 



For ease of calculation, we take 



i 

4>Um(x) = —=F k t{x)Y lm (Q) 

V 7TK 



u A =xi m {x c ) = -^/T(t)T{l + 2) S e (r c )Y em {Sl) (66) 
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The integral over dVt is immediate, giving 5 mm >^Wi and we suppress these delta functions in the following. Pulling 
out normalization factors, the inner product becomes 

r/ 2 _ _ 

aum = -A(£, k) / dt c cost c [Se(r c ) (d rc F k e{Ti,rf)) - (d rc Si(r c ))Fki{r),r f )] \ rc= o , (67) 

J-tt/2 



where A(£,k) = + 1)T(£)/V^k. 

For M = Q, Si(r c ) is independent of space t c and can be pulled out of the integral to give 

ri 2 

■ = -A(£,k){S e (r c )d rc - d rc S e (r c )} / dt c cos t c F M (rj, r/)| rc=0 

J-tt/2 



(68) 



Taking the limit r c — + (and remembering that Pii^" 1 has imaginary argument) one has ( [^5| , equation 8.1.4) 



S t (0) 



-3/2 

P^7o 1/2 (isinhr c ) , ,„ s o-i-i/2 /= 



-3/2 



yi cosh r c 
/2 B(£) 



r(^)r(^) 



(69) 



and 



^r c 5«(r c )| re= o = <9 r 



P_3y 2 1 ^ 2 (i sinh r c 



\/ 1 cosh r c 



-1-1/2, 

-1/2 



(i sinhr c ) 
(icoshr c ) 3 / 2 



r c =0 

e -3 l7 r/4 e 7r l (£/2+l/4)£ 2 _^-l/2 



(r(^)) 2 



o ™0>/2-1/2) 



(70) 



The inner product now becomes 



a Mm = -i l A(£,k){B(£)d rc +iC{£)} / dt c cost c Fki(r),r f )\r e =o 

-7T/2 



tt/2 



(71) 



where B(^), C(£) are all real. 

Again, both i«a=i £m and 0/t£ m are positive-frequency mode functions for the Bunch-Davies vacuum, so that 
(defined in equation vanishes: 



= (u A= ie m ,4> ke , ml ) 



= -i l A{l,k){-l) m+1 8 m -,-m{B{l)d rc + iC{£)} / dt c cast c F M (r,,r f )\ ra = 



ir/2 



The complex conjugate of this equation implies 

tt/2 



B(£)d ra dt c co8t c F ki (r},rf)\ rc=0 = iC(£) dt c cost c F ke (j],r f )\ rc=0 



tt/2 



-tt/2 



tt/2 



tt/2 



(72) 



(73) 



Defining 
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r /2 — 

Ii(k) = / dt c co8t c F k e(T),rf)\ rc =o (74) 



and using that je(kr) = (—l) je(—kr) is real, 



Thus we can write 



F/., (./■,■ I -/,))|,, ,,= ( 1 - i sin (_ f e < "°(-'°)j < (fccot(-t e )) 

1 + *— ^— )e^^fe)(-l)^(fccot(t c )) (75) 
sin(t c )/ 

= (-l)^*,^))!^ . 



9 rc / # c cost .F w (»7,r / )| re= o = i-=)^^(A) + (-l)^(fc)) . (76) 

-7T/2 



Substituting in, the full inner product is thus 



(%=ita,fe m ) = -2i' +1 C(*)A(*,fc)(Z*(*0 + (-l)%{k)) 



(77) 



2 

The calculation of the inner product thus requires the integral 



-{i t {k) + (-i)^(fc)). 



W 2 

It(k)= dtcCOBtcil + ibfie-^jiikrf)]^. (78) 







Changing coordinates and expressing 77, i c in terms of rf (there is only one free coordinate as r c has been fixed to 
zero), using 



1 



7? = - A /r| + 1 , tani c = -f-, (,!)') 



gives 



im = dr f {r ,l f 1)a/2 (l - ik^l) e ik V^ ji{k r f ) 

To proceed, note 



(80) 



ik / t*"^ I 1 

" (i-^ v / ^TT) e 4fe v^ T = -a r/ (i 7 l = ). (si) 



( ^ + i)3/ 2 v v / y~ V r / + 1 

For ^ odd, we need the imaginary part of Ii(k) and so can use ( (25|, equation 10.1.45) 
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1 Sill k \ I I r ~\~ .1. 

=- =iJ2^ri + l)j n (kr f )j n (k)P n (0) (82) 



kjrj + l 



Consequently, for I odd the integral of interest is 

00 />oo 

\mX lMd {k) =-ife^(2n + l)i„(A;)P n (0) / dr f d rf {j n {kr f ))j e (kr f ) (83) 



n=0 " u 

Using Mathematica, 

' . v \ £(£+l)-n(n + l) *-n 

dro r ] n (kr fjipikrf) = - — — — cos \n 

io j j\ £(f + 1) - n(n + 1) . ir£ . ?rn , Q .. 

(^ odd) = — --7- ry- f— sin — sin — (84) 

v ; (£-n-l)(£-n+l)(£ + n)(2 + £ + n) 2 2 y ' 

where it appears that n odd is required as well. However, looking at the sum (equation |8^), we see this term is 
multiplied by P n (0) which vanishes for n odd. So the only possibility for a nonzero term is if the denominator in this 
expression vanishes, that is if I = n ± 1. Substituting in these values and the definition of the Legendre polynomials 
we get 

—£ 7T 

Iml e , od d(k) = -ik-^—-P £ - 1 (0)(j^ 1 (k) +je+i(k)) 

= U^Pi-x(0)ji(k) (85) 

Including the prefactors (equation (|7§)) and after some algebra with T functions one gets that the overlap for £ odd 



OLktr, 



-i^/2£{£+l)k- 1 / 2 j l {k) ,£odd (86) 



agreeing with equation (p0[) . 

For £ even the calculation is more involved because the identity required in this case is ( [^5|, equation 10.1.46) 

cosfc./^ + l - Er=o(2" + l)yn(k)j n (kr f )P n (0) , 77 < 1 

7=— = (87) 

*V r / + 1 - EZo( 2n + l)y n (kr f )j n (k)P n (0) , r f > 1 . 

so that the range of integration in Tg(k) is split from < r f < 1 and 1 < rf < 00. As the resummation over akim^kim 
in equation ( filf ) works equally well for £ odd and even, it does not seem enlightening to pursue the calculation for 
even £ in all generality here. 

We can check a specific case, £ = 2, by integrating by parts and using the definition of ^(fcr/), 

2 . 



j 2 (kr f ) = (kr f Y --—d, 



1 \ sinfcr/ 



k 2 rf J kr f 



to get 
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yoo r r ( 1 — ik* /: 
J 2 (fc) = / dr, V V 



o 



lim fc 



(At/) 2 



1 ^ \ 

k 2 rj rf ) krf 



2 + i sin kr } 



(rj + l) 3 / 2 



r f 



k — 2" $t t 



r 2 f [1 - ik. rl + 1 , 

/V V f J ik^rj + 1 



(r 2 +1)3/2 



sinfcr/lg 



dr f 



sin /cr 



1. 



yjf f 1 L// f 



rj(l-ik^l) 
(rj + 1) 3 / 2 6 



+1 



(89) 



The only nonzero boundary term is the second one, at rj = 0, which can be read off, as it is only nonzero when the 
derivatives act on r/ rather than on ./rj + 1. Substituting r/ — > =F(a; — x _1 )/2 and using Mathematica gives 



(90) 



(91) 



„. 6i(e lk -l) 6e ik (1 + 2e lk ) 



Combining the integral with the prefactors (eqn. (78)), and taking the real part gives 

[3 1 

ctk,i=2.m — — z2 W — -^-g [3 sin k — 3k cos k — k 2 sin k\ = —2i\[3 kT 1 ! 2 ji(k), 
again in agreement with equation j5C|). 

C. Overlap along r\ = 0, for I — 1 



Here we discuss the integral (equation (|3 

7/(*)|E 



*7 [ r / 



evaluated on the surface in region F corresponding to 77 = 0. For 77 = 0, g = (1 + rj)/(2rf) and so 



= / dr f r f 2 j e (kr f )(g 2 - 1)0^) 
Jo 



becomes for £ = 1 

h(k)\ v=0 = 
This can be rewritten as 
h(k)\ v=Q = 



2 1 



tt(3) y 



sinfcrj cos kr 



(kr f ) 2 krf 
sm krf cos krf 



8rr 



. 2 



(i-^r ^ Ai+^r A 1+*-/ 



2ry 



(krf) 2 krf 



(l-?7 2 ) 



2N2 



In 



1 + r/ 



l-»7 



1 + ry 



(92) 



(93) 



(94) 



(95) 



and integrated using Mathematica by considering 
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+ 



f dr 
Jo 

/>OC 

I drf 



sin At/ cos At/ 



(At/) 2 
sin At/ 



At/ 
cos At/ 
krt 



In 



1-r/ 

l + ?7 
-1 + r/ 



1 + 



1 + r/ 2 



(96) 

(97) 
(98) 



to get 



This gives 



/ 1 4 2 

7i(fc) = -\ — { lim — — sin At/ + — [-(ci(oo) - ci(-oo))(cos(fc) + A; sin fc) 

V 27T r/^oo 3fc z k° 

2lTT 

+ (Si(oo) — Si(— oo))(— A; cos Ac + sin Ac)] + -^-(cosA; + A; sin Ac)} 

1~ / 4 2ir \ 

— I lim — -rsinArf-l — — =- (sin fc — fc cos fc) ) . 
2tt Vr/-oo 3fc 2 J fc 5 / 



Q!fe,<=l,m = 1 



( lim — ^-sinAr/H — ^ (sin fc — fc cos Ac)) = — 2ik 1 / 2 ji(fc), 

yr f ^oo 3fc^ fc b / 



(99) 



(100) 



where the oscillating first term has been dropped as it has the wrong asymptotics as fc — > oo. (If lim r/ ^ 00 sinfcr/ was 
finite rather than the giving zero in this limit, the overlap would diverge as fc — > oo.) With this, aum again agrees 
with the results found above. 



ACKNOWLEDGMENTS 



D.K. thanks A. Guth, and J.D.C. thanks A. Anderson, S. Axelrod, L. Ford, D.E. Freed, A. Kent, and A. Vilcnkin for 
conversations and is especially grateful to M. White. J.D.C. also thanks the Aspen Center for Physics, A. Loeb of the 
Harvard-Smithsonian Center for Astrophyics, the Tufts Institute of Cosmology and the Insitute d'Astrophysique de 
Paris for hospitality in the course of this work. The work of D. K. has been supported in part by NSF PHY-92-18167. 
J. D. C. is supported by an NSF Career Advancement Award, NSF PHY-9722787, and, at the commencement of this 
work, by an ONR grant as a Mary Ingraham Bunting Institute Science Scholar at Harvard University. 



[1] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge University Press, New York, 1982). 

[2] S. A. Fulling, Aspects of Quantum Field Theory in Curved Spacetime (Cambridge University Press, New York, 1989). 

[3] A. H. Guth, Phys. Rev. D 23, 347 (1981). See also E. W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, 

New York, 1990); A. D. Linde, Particle Physics and Inflationary Cosmology (Harwood, Chur, 1990). 
[4] M. Bucher, A. Goldhaber, and N. Turok, Phys. Rev. D 52, 3314 (1995); Nucl. Phys. Proc. Suppl. 43 (1995) 173 . 
[5] K. Yamamoto, M. Sasaki, and T. Tanaka, Astrophys. J. 455, 412 (1993); A. D. Linde, Phys. Lett. B351, 99 (1995); A. 

D. Linde and A. Mezhlumian, Phys. Rev. D 52, 6789 (1995). 



21 



S. Coleman and F. de Luccia, Phys. Rev. D 21, 3305 (1980); J. R. Gott, Nature 295, 304 (1982); J. R. Gott and T. S. 
Statler, Phys. Lett. B136, 157 (1984). 

D. H. Lyth and A. Woszczyna, Phys. Rev. D 52, 3338 (1995); J. Garcia-Bellido, A. R. Liddle, D. H. Lyth, and D. Wands, 
Phys. Rev. D 52, 6750 (1995). 

M. Sasaki, T. Tanaka, and K. Yamamoto, Phys. Rev. D 51, 2979 (1995). 



U. Moschella and R. Schaeffer, "Quantum Fluctuations in an Open Universe," preprint gr-qc/9707007 . 

A. D. Linde, Phys. Lett. B116, 335 (1982); A. A. Starobinsky, Phys. Lett. B117, 175 (1982); A Vilenkin and L. H. Ford, 

Phys. Rev. D 26, 1231 (1982); B. Allen, Phys. Rev. D 32, 3136 (1985); L. H. Ford and A. Vilenkin, Phys. Rev. D 33, 2833 

(1986); B. Allen and A. Folacci, Phys. Rev. D 35, 3771 (1987); D. Polarski, Phys. Rev. D 43, 1892 (1991). 

K. Kirsten and J. Garriga, Phys. Rev. D 48, 567 (1993). 

A. D. Linde and A. Mezhlumian in reference 5 above, J. Garcia-Bellido, Phys. Rev. D 54, 2473 (1996) 
J. Garriga, Phys. Rev. D 54, 4764 (1996). 

K. Yamamoto, M. Sasaki, and T. Tanaka, Phys. Rev. 54, 5031 (19 96); M. Sasaki, T . Tanaka, Phys.Rev. D 54 4705 (1996). 



J. Garcia-Bellido, J. Garriga, X. Montes, "Quasi-Open Inflation," hep-ph/9711214 
D. I. Kaiser, in preparation. 

L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. Lett. 73, 3195 (1994); Y. Shtanov, J. Traschen, and R. Branden- 
berger, Phys. Rev. D 51, 5438 (1995); D. Boyanovsky et al, Phys. Rev. D 51, 4419 (1995), 52, 6805 (1995), and 54, 7570 
(1996); M. Yoshimura, Prog. Theo. Phys. 94, 8873 (1995); D. I. Kaiser, Phys. Rev. D 53, 1776 (1996). 
J. D. Cohn and D. I. Kaiser, in preparation. 

V. Mukhanov, H. Feldman, and R. Brandenberger, Phys. Rep. 215, 203 (1992); V. Mukhanov, L. Abramo, and R. 
Brandenberger, Phys. Rev. Lett. 78, 1624 (1997) and Ph ys. Rev. D 56 , 3248 (1997); See also M. Sasaki, T. Tanaka, 



"Super-Horizon Scale Dynamics of Multi-Scalar Inflation," gr-qc/9801017 
B. Allen, Phys. Rev. D 51, 3136 (1985). 

S. W. Hawking and G. F. R. Ellis, Large-Scale Structure of Spacetime (Cambridge University Press, New York, 1973). 
T. S. Bunch and P. C. W. Davies, Proc. Roy. Soc. A 360, 117 (1978); T. S. Bunch and P. C. W. Davies, J. Phys. A 11, 
1315 (1978). See also R. H. Brandenberger, Nucl. Phys. B245, 328 (1984); A. H. Guth and S.-Y. Pi, Phys. Rev. D 32, 
1899 (1985). 

M. Bucher and N. Turok, Phys. Rev. D 52, 5538 (1995). 

B. S. DeWitt, in Relativity, Groups, and Topology II, Les Houches 1983, ed. B. S. DeWitt and R. Stora (North-Holland, 
New York, 1984). 

M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (Dover, New York, 1965). 

See for example T. Hamazaki, M. Sasaki, T. Tanaka, K. Yamamoto, Phys. Rev. D 53, 2045 (1996); J. Garcia-Bellido, 
Phys. Rev. D 54, 2473 (1996); J. Garriga, Phys. Rev. D 54 4764 (1996); J. Garriga, X. M ontes, M. Sasaki, T. Tanaka, 



"Canonical Quantization of Cosmological Perturbations in the One Bubble Open Universe," istro-ph/970622£ . 
R. Rajaraman, Solitons and Instantons (North-Holland, New York, 1982). 

I. S. Gradshteyn, I. M. Ryzhik, Table of Integrals, Series and Products, Fourth Edition (Academic Press, Inc., San Diego, 
1980). 



22 



